How a liquid becomes a glass both on cooling and on heating 
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The onset of structural arrest and glass formation in a concentrated suspension of silica nanoparticles in a 
water-lutidine binary mixture near its consolute point is studied by exploiting the near-critical fluid degrees of 
freedom to control the strength of an attraction between particles and multispeckle x-ray photon correlation 
spectroscopy to determine the particles' collective dynamics. This model system undergoes a glass transition 
both on cooling and on heating, and the intermediate liquid realizes unusual logarithmic relaxations. How 
vitrification occurs for the two different glass transitions is characterized in detail and comparisons are drawn to 
theoretical predictions for glass formation in systems with attractive interactions. 
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The glass transition remains a grand challenge for con- 
densed matter science [1 J2|], with applications in diverse areas 
such as protein folding [3|] and the flow of granular materials 
lit]. Especially remarkable then are predictions Q 0, 0, 01, 
that in systems with short-ranged attractions a dense liquid 
may become a glass both on cooling and on heating. The sce- 
nario is as follows: When the ratio of attraction strength to 
temperature is small, the attractions play no role and, at hig h 
densities, a repulsion-dominated glass (RG) is realized fl fioll . 
However, as the ratio is increased, surprisingly, attractions 
cause the glass to melt. With still further increase, the resul- 
tant liquid re-vitrifies into an attraction-dominated glass (AG). 
The intermediate liquid is itself highly unusual in that its den- 
sity fluctuations are predicted to relax logarithmically in time. 
In the present paper, by exploiting fluid-mediated interactions 
between silica particles in a near-critical water-lutidine binary 
mixture to control the attraction strength 01 U Il2l 1 1 311 - and 
multispeckle x-ray photon correlation spectroscopy (XPCS) 
lfl4l [I5I [l6ll to determine the collective dynamics, we show 
that re-entrant glassy behavior and logarithmic relaxations are 
realized in this model system, simply by changing tempera- 
ture. We are thus able to characterize in unprecedented detail 
how the non-ergodicity parameter of the RG decreases upon 
attraction-driven melting, and how relaxations in the interme- 
diate fluid arrest as the attactive glass phase is approached. 

Glassy states in colloidal suspensions are essentially spa- 
tially and temporally inhomogeneous |l7 , 18 ] . It is therefore 
paradoxical that mode-coupling theory (MCT) - a mean-field- 
like approach - is the only theory to-date capable of quantita- 
tive predictions concerning collective dynamics near the glass 
transition ifToh . Nevertheless, MCT captures essential aspects 
of the dynamical arrest observed near the prototypical glas s 
transition that occurs in a suspension of hard spheres J9I I10I1 . 
where arrest occurs at sufficiently high densities, when each 
particle becomes trapped in a cage of its neighbours. 

Instead of caging, it is intuitive that sufficiently strong in- 
terparticle attractions permit an alternate route to dynamical 
arrest, in which particles simply stick together. Recent MCT 
calculations earned out to investigate this possibility lead to 
the novel predictions of Refs. Support for these pre- 

dictions comes from optical PCS experiments 
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in which a depletion attraction between particles was created 
by the addition of non-adsorbing polymer. However, how 
the behavior evolves with attraction strength was only deter- 
mined with limited resolution, in part because tuning the at- 
traction requires laboriously mixing samples. By contrast, sil- 
ica particles in near-critical water-lutadine (WL) mixtures ex- 
perience a temperature-dependent attraction as a result of the 
confinement of the fluid critical fluctuations between particles 
01211 . and evidenced by the particles' reversible aggregation 
as the temperature is increased towards the two-phase region 
01 ill . Recently, Pontoni et al. ill 311 confirmed the temperature- 
dependent attraction via small-angle x-ray scattering (SAXS) 
and XPCS in a dilute sample. The present study was carried 
out to explore in detail how glassy behavior depends on attrac- 
tion strength in concentrated suspensions of silica particles in 
near-critical WL mixtures. 

Illustrated in Fig. Q2a) is the SAXS from 195 nm-radius 
silica spheres (polydispersity < 10%) at a volume fraction of 
<j) ~ 0.52 in WL with a lutadine concentration of c = 0.24. Par- 
tially coherent illumination gives rise to a speckled intensity, 
whose mean and time-autocorrelation characterize the sam- 
ple's static structure and dynamic relaxations, respectively. To 
this end, sequential SAXS patterns were collected using CCD- 
based area detectors at frame rates of up to 330 Hz (Run A and 
Run B) or 60 Hz (Run C) C3]. Fig. [flb) shows the time and 
azimuthally averaged SAXS intensity (symbols) plotted ver- 
sus wavevector (Q) in the neighbourhood of the peak in the 
interparticle structure factor, obtained at several temperatures 
during Run A. At the lowest temperature shown (33.400°C), 
there is a relatively narrow peak at QR = 3.6. With increasing 
temperature, the peak gradually becomes broader and shifts 
slightly toward larger QR, indicating a decrease in the cor- 
relation length, accompanied by slight decrease in the mean 
particle spacing. Concommitantly, the intensity increases for 
QR < 2.5, corresponding to the emergence of larger scale 
inhomogeneities. To quantify this evolution, we have fitted 
these data to the structure factor of spheres subject to a po- 
tential consisting of a hard core repulsion together with an at- 
tractive square well within the mean-spherical approximation 
|0]. The best fits - the solid lines in the figure - satisfactorily 
reproduce the main features of the data, especially in view of 
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FIG. 1: (a) False-color image of the SAXS pattern from concentrated 
silica spheres in WL, obtained at beamline 8-ID at the Advanced 
Photon Source (APS) in Run A with a 0.15 s exposure, (b) Average 
SAXS intensity vs. wave vector obtained for several temperatures on 
warming (Run A). The lines are the best fits to the model described 
in the text. 



the fact that the only parameters varied in each fit were the 
overall intensity and the strength of the attraction (U). The 
particle radius and the hard core radius (R = 200 nm) were de- 
termined from fits to data obtained at the lowest temperature, 
where the attractive component is negligible. Because the 
model lineshape depends only on the dimensionless combina- 
tion u = UA/(IcbTR) within the g-range studied, the width (A) 
of the square-well attraction was fixed to be 0.017?. At each 
temperature, the silica volume fraction was determined via 
the corresponding x-ray transmission of the sample. Fig. [3a) 
shows the best fit values of u for run A (diamonds), and two 
additional warming runs: B (triangles) and C (circles). In each 
case, u takes on a small value at the lowest temperatures, but 
increases rapidly at a certain onset temperature. We ascribe 
the different onset temperatures for the different runs to small 
variations in the lutidine concentration between runs. 

Normalized intermediate scattering functions (g\) 11611 at 
QR = 2.34 for delay times from 17 ms to 410 s for a sam- 
ple with (]) ~ 0.49 (Run C) are illustrated for temperatures 
from 33.000°C to 33.375°C in Fig.EJa), and from 33.390°C 
to 33.600°C in Fig. [2b). At the lowest temperature, g\ shows 
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FIG. 2: (a) Best fit values of u for Run A (diamonds), B (triangles) 
and C (circles). The different colors indicate the different phases 
realized: red corresponds to a RG, blue to a liquid, and green to 
an AG. The fitted value of u appears constant in the AG, suggest- 
ing that no structural evolution occurs in this phase, (b) Phase di- 
agram in the u — ty plane. Diamonds, triangles, and circles corre- 
spond to data points from Runs A, B, and C, respectively. Red, 
green, and blue symbols correspond to RG, AG, and liquid phases, 
respectively. Solid lines are schematic phase boundaries, inspired by 
Refs. |5j, |a |7|, la, with the black circle corresponding to the MCT 
A3 singularity. The dashed line separates the region of logarithmic 
relaxations above the line, from the region of stretched exponential 
relaxations below. Although our particles are subject to gravitational 
settlement, the volume fraction across the height of the x-ray beam 
(20 ,um) was approximately constant to within 1%. Usually, the vol- 
ume fraction also remained constant in time. In Run B, however, the 
volume fraction varied as a result of slow settlement. 



an initial rapid decay (which is not fully captured by the 60 Hz 
frame rate of the camera employed in Run C) to reach a value 
of about 0.5 at 0.1 s. For times between 0.1 s and 100 s, g\ re- 
mains nearly constant, revealing arrested dynamics over three 
decades in time. Such behavior is similar to that observed for 
hard-sphere glasses and is usually interpreted in terms 

of caging. The plateau value (/) is called the "non-ergodicity 
parameter". Because of the limited time range probed in the 
present study (< 410 s), we cannot determine whether the de- 
crease in g\ observed for times longer than about 200 s is a 
long-time equilibrium relaxation or is associated with aging 
of an arrested state l23ll . As the temperature is increased, / 
decreases continuously to reach zero at 33.375°C [Fig. O a)], 
while the initial decay continues to about 0.3 s. Thus, within 
this range of tempertures, for which the interparticle attrac- 
tions are increasing [Fig. |2fa)], remarkably, the initially- 
arrested RG melts into a dynamic liquid. The observed evo- 
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FIG. 3: gi versus delay time during Run C. For clarity, (a) cor- 
responds to temperatures from 33.000°C to 33.375°C and (b) from 
33.390°C to 33.600°C. Points are data. Lines correspond to the mod- 
els discussed in the text. 
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lution of g\ is distinct from what occurs when a hard sphere 
suspension undergoes the liquid-glass transition as a function 
of volume fraction. In that case, a plateau with a non-zero, 
Q-dependent value starts to emerge within the liquid phase, 
and, with further increase in the volume fraction, the plateau 
extends to longer times, showing only a small increase in its 
value through the glass transition j^l fioll . 

On warming beyond 33.375°C [Fig. Ob)], the relaxation 
of gi becomes progressively slower with increasing temper- 
ature. On reaching 33.485°C, for times from about 0.1 s to 
about 100 seconds, g\ appears nearly linear in the semi-log 
plot of Fig. Ob), indicating the onset of nearly logarithmic 
relaxations. On warming still further, the logarithmic behav- 
ior continues, but with a slope that decreases with increasing 
temperature, so that by 33.600°C, g\ barely decreases to 0.95 
after 300 seconds, corresponding to an almost completely ar- 
rested state. As noted previously, logarithmic behavior of gi 
is expected for the liquid near an AG. 

To quantify these observations, for temperatures 
from 33.000°C to 33.480°C, we have fitted gi us- 
ing an empirical, modified-stretched-exponential form: 



/+(!-/)« 



-C/t ) a 



To limit the num- 
ber of parameters varied in the fits, the stretching ex- 
ponent (a) was fixed to 0.6, leaving /, to and ti as 
fitting parameters. From 33.485°C to 33.600°C, we 
have instead modeled g\ via an alternate empirical form: 
gl = [e-'/"+(l-e-'/")[l-51og(l+r/fi)]]e-('/'2) 3/2 , 
which manifests logarithmic relaxations for t% ^> t 3> t \ . On 
the semi-log axis of Fig. Ob), S is the slope of the logarithmic 



FIG. 4: (a) /, (b) to, and (c) 5 vs. temperature for Run B at QR = 2.28 
(cyan squares), 2.86 (blue diamonds) and 3.70 (orange triangles), 
and for Run C at QR = 1.13 (red circles), 2.34 (magenta squares), 
and 2.82 (green diamonds). The solid lines are guides-to-the-eye. 



gl. In this case, the fitting parameters are S and t%, while the 
onset of the logarithmic behavior (fi) was fixed at 0.02 s. 
The solid lines in Figs. O a ) and Ob) correspond to the best 
fits to these models, providing a satisfactory account of the 
measured g\ at all temperatures, except at 33.480°C, where 
gl is evolving from one form to the other. In both equations, 
the compressed exponential factor Il23ll . involving t2, accounts 
for the behavior at long times (?2 > 200 s) and will not be 
discussed further in this paper. 

Fig. Ufa) shows the best-fit values of / plotted versus tem- 
perature, obtained during run B at QR = 2.28, 2.86, and 3.70 
and during Run C at QR = 1.13, 2.34, and 2.82. As may be 
seen, / takes a (^-dependent value with systematically larger 
values near the peak of the static structure factor (QR ~ 3.7). 
This is similar to what is observed for hard-sphere glasses. 
However, here, / goes continuously to zero with increasing 
temperature. This behavior stands in contrast to what oc- 
curs at the melting of a hard-sphere glass with decreasing 
density, where the non-ergodicity parameter goes to zero dis- 
continuously. MCT also predicts a discontinuous decrease in 
the non-ergodicity at the RG-to-hquid transition in the pres- 
ence of attractive interactions |5l |6j lTJ. Thus, this aspect of 
our measurements presents a contradiction to MCT. It may 
be noted, though, that certain MCT models [24] do permit a 
glass transition involving continuous increase in / as we ob- 
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serve in Run C. For Run B, density correlations relax very 
slowly, so that g\ does not decrease to zero for times less than 
200 s at any temperature. In fact, g\ relaxes progressively 
more rapidly for temperatures increasing to 33.55°C, but re- 
laxes progressively more slowly for temperatures increasing 
beyond 33.55°C. This leads to an apparent minimum in the 
fitted value of / at 33.55°C, before g\ crosses over to the loga- 
rithmic form for 33.64°C and above. We ascribe the difference 
in behavior between Run B and Run C to the closer proximity 
of a glass-to-glass transition in Run B. 

The best-fit values of the relaxation time (fo) are shown in 
Fig. |4fb) plotted on a logarithmic scale versus temperature. 
Evidently, fo depends strongly on temperature. For both Runs 
B and C, fo increases rapidly with temperature from on the 
order of 0.02 s to reach about 0.5 s, before the relaxations 
become logarithmic, fo depends relatively weakly on Q, al- 
though in Run B it is somewhat larger at the structure fac- 
tor peak (QR = 3.7) than elsewhere. For Run C, it was not 
possible to determine fo at the structure factor peak because 
of the presence of Bragg peaks from small coexisting crys- 
tallites. In comparison, a Q-independent time-scale, that in- 
creases rapidly with density, is a hallmark of relaxations near 
the glass transition of hard spheres. 

The best-fit values of S within the logarithmic regime are 
summarized in Fig. Uc) versus temperature. S is only weakly 
g-dependent. Interestingly, near the temperature at which it 
vanishes, S appears to decrease approximately linearly to zero 
with increasing temperature. Although expressions for S are 
available within the context of MCT for small displacements 
in u — (]) phase space from the so-called A3 singularity point 
1 2411 . it is unclear how these predictions relate to the behavior 
seen in Fig.Uc). 

Our results for the phase behavior of silica particles in WL 
in the volume fraction-interaction strength ((j) — u) plane are 
summarized in Fig. Efb), where diamonds correspond to Run 
A, triangles to Run B, and circles to Run C. Each point is plot- 
ted with u determined by fitting the SAXS lineshape and with 
(j) determined by the corresponding x-ray transmission. Red 
symbols correspond to the RG, defined for the purposes of the 
figure by the criterion that g\ > 0.17 at t = 10 s. Green sym- 
bols correspond to the AG, defined similarly via g\ > 0.81 at 
f = 10 s. Blue symbols correspond to the liquid. In the con- 
text of Refs.jsllaSH it is natural to add the schematic phase 
boundaries shown as solid lines in the figure. The dashed line 
in the figure separates the region of logarithmic relaxations 



above the line from the region of stretched exponential relax- 
ations below. Overall, there is good agreement between our 
experimental phase diagram and that predicted. In addition, 
this study reveals a number of previously unremarked aspects 
including: (1) A continuous decrease in the non-ergodicity 
parameter with attraction strength through attraction-driven 
melting, (2) the existence of a weakly ^-dependent time-scale 
that increases through attraction-driven melting, and (3) the 
continuous decrease of S for temperatures increasing towards 
the liquid-to-attractive glass transition. We thank T. Chiba, 
E. R. Dufresne, R. L. Leheny, C. O'Hern, S. Sanis, M. Span- 
nuth, M. Trias, and J. Wettlaufer for discussions, a referee for 
key clarifications, and the NSF for support via DMR 0453856. 
The APS is supported by the DOE. 
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